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X-rays from Colliding Stellar Winds: the case of close 
WR+O binary systems 
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ABSTRACT 

We have analysed the X-ray emission from a sample of close WR+O binaries using 
data from the public Chandra and XMM-Newton archives. Global spectral fits show 
that two-temperature plasma is needed to match the X-ray emission from these objects 
as the hot component (kT > 2 keV) is an important ingredient of the spectral models. 
In close WR+O binaries, X-rays likely originate in colliding stellar wind (CSW) shocks 
driven by the massive winds of the binary components. CSW shocks in these objects 
are expected to be radiative due to the high density of the plasma in the interaction 
region. Opposite to this, our analysis shows that the CSW shocks in the sample of 
close WR+O binaries are adiabatic. This is possible only if the mass-loss rates of the 
stellar components in the binary are at least one order of magnitude smaller than 
the values currently accepted. The most likely explanation for the X-ray properties of 
close WR+O binaries could be that their winds are two-component flows. The more 
massive component (dense clumps) play role for the optical/UV emission from these 
objects, while the smooth rarefied component is a key factor for their X-ray emission. 

Key words: shock waves - stars: individual: WR46, WR47, WR79, WR139, WR141, 
WR145, WR148 - stars: Wolf-Rayet - X-rays: stars 



1 INTRODUCTION 

As first proposed by IPrilutskii fc Usovl |l976) and 
ICherepashchuk I (1 19761 ). the present concept on the origin 
of X-rays in massive Wolf-Rayet+O (WR+O) binaries as- 
sumes that they arise in colliding-stellar-wind (CSW) shocks 
resulting from the interaction of the massive winds of the 
stars in the binary. High wind velocities and high mass-loss 
rates of WR and O stars suggest that the CSW region will 
be a luminous X-ray source. Thus, WR+O binaries are ex- 
pected to be brighter (more luminous) in X-rays than single 
WR stars an d, in fact, thi s was found in the early surveys 
of WR stars (|PollockUl987h . 

It is interesting to note that wide WR+O binaries are 
among the brightest WR objects in X-rays and their grat- 
ing spectra with the modern observatories ( Chandra, XMM- 



dicate thermal emission (Skinner ct al. 2001 


; Raasscn ct al. 


2003; 


Schild et al.ll2004|Pollock et all 20051; 


Zhekov & Park 


2010b 


). Analyses of these data confirmed that X-rays from 



CSW shocks is the most probable physical picture for the X- 
ray emission from wide WR+O binaries (although in some 
cases the observed pict ure might be more complex than that; 
IZhekov fc Parkll2010al ,bb 
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On the other han d, it is noted from those early X-ray 
surveys (|Pollocklll987h that close WR+O binaries are not 
as luminous in X-rays as the wide WR+O binaries. This 
is an important observation since the stellar wind param- 
eters (which in general determine the X-ray emission from 
CSWs) are practically within the same range for both types 
of binary systems. 

As a possible resolution to this problem. ICherepashchukl 

proposed that most of the mass flux in the stel- 
lar wind of a WR star is in the form of dense clouds. 
These clouds penetrate freely the CSW region formed by 
the continuous component of the winds and they do not 
contribute to the X-ray emission from the binary. Thus, in a 
close WR+O binary only small part of the stellar wind gas 
(the continuous component) plays role for the X-ray emis- 
sion from the system. However, in a wide WR+O binary 
the CSW region forms far from the stars where the dense 
clouds have already expanded and 'dissolved' in the contin- 
uous component of the wind, thus, the entire mass flux of 
the stellar wind contributes to the X-ray emission from the 
binary. 

Since the modern X-ray observatories ( Chandra, XMM- 
Newton) provide us with data of much higher quality (both 
in sensitivity and spatial resolution) compared to that of 
previous X-ray missions, we decided to explore the case of 
CSWs in close WR+O binaries in some detail. This was 
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done for a number of such objects that have been observed 
by Chandra and/or XMM-Newton. We define the sample of 
close WR+O binaries in Section [2] and we describe the cor- 
responding X-ray data in Section [3] In Section [4] we present 
the results from the global spectral models. In Section [5] we 
discuss our results and we list our conclusions in Section [5] 



2 THE CSW BINARIES SAMPLE 

In the sample of close CSW bin aries, we included th ose ob- 
jects of known WR+O binaries l|van der Huchtll200ll ; Tables 
18 and 19) with orbital period smaller than 20-22 days. We 
have checked the archives of the modern X-ray observato- 
ries {Chandra and XMM-Newton) for available data on the 
thus defined sample. Some objects with X-ray detection are 
not included in our study and we particularly mention the 
following objects: WR48 (because it is a suggested triple 
stellar system and it s X-ray emission is quite complex; see 
ISugawara et alfeOOST ); WR43a (since it is in the core of NGC 
3603, thus, the object cannot be resolved); WRIOlk (it is in 
the center of the Galaxy and unresolved: being in l'.'8 from 
Sgr A; SIMBAD). 

Table 1 presents the list of close CSW binaries consid- 
ered in this work and summarizes some of their properties. 



3 X-RAY DATA 

Since the studied objects are not very bright in X-rays, we 
made use of the available data from the Advanced Cam- 
era for Imaging Spectroscopy (ACIS-S) on-board Chandra 
and of the data from the European Photon Imaging Cam- 
era (EPIC; pn and MOS detectors) on-board XMM-Newton. 

For the analysis of the ACIS-S data, we extracted the 
X-ray spectra following procedures in the Science Threads 
for Imaging Spectroscopy of the CIAO Q 4.3 data analy- 
sis software. The response functions and ancillary response 
functions for all spectra were generated using the Chandra 
calibration data base CALDB v4.4.2. 

For extracting X-ray spectra from the EPIC (pn, MOS) 
data, we made use of the XMM-Newton SAS 10.0.0 data 
analysis software. The SAS pipeline processing scripts em- 
proc and epproc were executed to incorporate the most re- 
cent calibration files. The data were then filtered for high 
X-ray background following the instructions in the SAS doc- 
umentation (e.g., adopting typical threshold rates for high- 
energy background in pn and MOS: 0.4 and 0.2 cts s~\ 
respectively). The SAS procedures rmfgen and arfgen were 
adopted to generate the corresponding response matrix files 
and ancillary response files for each spectrum. For each data 
set, the MOS spectra in this analysis are the sum of the spec- 
tra from the two MOS detectors. 

Table 2 presents the basic information about the X- 
ray data of the objects in our sample of close CSW binaries. 
Also, some details on the individual objects are given below. 

WR46. There is one pointed observation of this ob- 
ject with XMM-Newton. The EPIC data are of good quality 

1 Chandra Interactive Analysis of Observations (CIAO), 
|http: //cxc.harvard.edu/ciao/| 

Science Analysis Software (SAS),[http://xmm.esac.esa.int/sas/| 



(no long background flares), thus, all the spectra, pn and 
combined MOS, were included in this analysis. This is the 
object with the shortest binary period in our sample (see 
Table [TJ, thus, the EPIC spectra are average (integrated) 
over the binary period (also, see Appendix [X] for discussion 
of the variability issue). 

WR47. There are three pointed observations of this 
star in the XMM-Newton archive and only two of them (see 
Table r2J) were used in this study since the third one showed 
considerable high-energy X-ray background rates. Both pn 
and MOS spectra were included in the current analysis. 

WR139. There are six pointed observations of WR139 
carried on with XMM-Newton, and here we used the pn 
spectra of three of them that have good photon statistics. 
The basic reason for the lower quality of the rest of the 
data sets is the high X-ray background which shortens their 
effective exposure. 

WR141. There are two data sets in the XMM-Newton 
archive that detect this WR binary at ~ 5' 1 from the optical 
axis. Unfortunately, due to the high X-ray background rate, 
only part of the MOS exposure in the second observation 
could be used in our analysis. 

WR145. We find one Chandra and three XMM- 
Newton archive observations that have detected this WR 
star respectively at ~ 5!2 and ~ 9.'6 from the optical axis. 
The quality of the EPIC data is not high and only the MOS 
spectra could be used in this analysis. 

WR148. In the XMM-Newton archive, there are two 
pointed observations of this distant WR binary. Due to high 
X-ray background, the pn spectrum of the second data set 
was excluded from our analysis. 

WR79. We find two Chanrda and eight XMM-Newton 
archive observations respectively at ~ 1' and ~ 3' from the 
optical axis, but only the former provide useful X-ray spectra 
of this WR object. The reason is that there are two near- 
by sources (within 6" from WR79) that cannot be spatially 
resolved by XMM-Newton (see Appendix [A] for details). 



4 GLOBAL FITS 

Since it is expected that the X-ray emission from the close 
WR+O binaries originates in the interaction region of their 
massive winds, a systematic approach to modelling these 
spectra could be to use a hydrodynamic CSW model. How- 
ever, the quality of the current data (undispersed CCD spec- 
tra with not very good photon statistics for most of the spec- 
tra) is not a good basis for carrying on such a modelling in 
detail. This is why, we decided to use conventional, discrete- 
temperature plasma, models for the global fits to the X-ray 
spectra of the objects in our sample. Thus, we adopted the 
optically thin plasma model apec i n a recent vers ion (11.3.2) 
of the XSPEC analysis package (|Arnaudlll996i ). Since the 
shocked WR gas dominates the X-ray emission from CSWs 
in WR+O binaries (the mass loss of a WR star is about 
an order of magnitude higher than of an O star), in the 
global spectral fits we kept the abundances of the chemical 
elements fixed to their values typical for the WN and WC 
stars, correspondingly (|van der Hucht et al.lll986h . To de- 
rive the X-ray absorption towa rds each object, we used the 
iMorrison fc McCammor] l|l983h cross-sections (model wabs 
in XSPEC). And it is worth noting that we found X-ray 
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Table 1. Properties of the close WR+O binaries 



Name 


Spectral 


Ad 


Distance 


Period 




M 


Sources 




type 


(mag) 


(kpc) 


(days) 


( km s- 1 ) 


( M Q yr- 1 ) 


used for M 


WR46 


WN3+OB 


1.05 


4.07 


0.329 


2450 


9.94e-06 


2, 3 


WR47 


WN6+05 


3.96 


3.80 


6.24 


1660 


2.30e-05 


6, 7, 10, 12 


WR139 


WN5+0 


2.83 


1.90 


4.21 


1600 


1.70e-05 


1, 4, 5, 8, 10, 11, 12, 13 


WR141 


WN5+05 


4.10 


1.26 


21.7 


1550 


2.48e-05 


1, 8, 10 


WR145 


WN7/WCE+? 




1.70 


22.55 


1390 


3.68e-05 


1, 10 


WR148 


WN8+B3 


2.58 


8.28 


4.32 


1500 


3.95e-05 


8, 9, 12 


WR79 


WC7+05-8 


1.54 


1.99 


8.89 


2270 


5.59e-05 


1, 6, 10, 11, 12 



Note — The object name, spe ctral type, interstella r extinction, distance, orbital period and wind velocity are from the VH-th catalogue 
of galactic Wolf-Rayet stars llvan der H ucht 200 lj; A v = I.I IAy) with the following exceptions: (a) the win d velocity for WR47 is the 
mean value f or WN6 stars from Eenens &: Williams! l ll994h an d for WR148 is from [ Nishimaki et all l l2008h ; (b) the orbital period for 
WR46 is from lMarchenko et al.l [|2000l l and for WR145 is from lMuntean et alj l l2009t) ; fcl the distance to WR145 is that to Cyg OB2. 
The mass-loss rates are the mean values for each object as found in the literature with no clumping taken into account. The individual 
data were scaled correspondingly if derived at value s for the distance to the object and /or its stella r wi nd velocity different from thos e 
ad opted here. The c orres ponding references are: (11 lAbbott et ah| l|l986|l: (21 [ Crowther et alj lll 995h: ( 3l lHaman n fc Koesterke ] lll99gfk 
f4l lHirv et alj [|2006|1 ; (51 iKurosawa et alj J2002h; f 6 ULamontagne et alj 1U996T I (7 llMoffat et alj dl990fl: f 81 iNishimaki et alj <2008l) ; (9) 
iNugis et al.1 .] 199,4): flOl lNugis fc Lamerd [|2000l l; Qll lPrinia et alj jl99of >: (12) ISt- Louis et al.l l|l988h : (131 ISt- Louis et a.1.1 (fl99a ) . 

Table 2. X-ray Observations 



Name 


Observation 
IDs 


Binary 
phase 


Detector 


Source 
counts 


Rate 

( cts s _1 ) 


Data 


WR46 


0109110101 


0.55 


pn 


3790 


6.20e-02 


3, 4, 5 








MOS 


2605 


3.62e-02 




WR47 


0109480101 


0.62 


pn 


1911 


4.50e-02 


1 








MOS 


1665 


3.22e-02 






0109480401 


0.89 


pn 


1004 


2.86e-02 










MOS 


843 


1.99e-02 




WR139 


0206240201 


0.47 


pn 


1514 


0.176 


1, 2 




0206240701 


0.55 


pn 


1723 


0.214 






0206240801 


0.81 


pn 


3744 


0.273 




WR141 


0404540201 


0.02 


MOS 


734 


5.30e-02 


5 


WR145 


10969 


0.33 


ACIS-S 


757 


2.64e-02 






0165360101 


0.00 


MOS 


297 


2.06e-02 






0165360201 


0.09 


MOS 


366 


2.54e-02 






0165360401 


0.18 


MOS 


361 


2.50e-02 




WR148 


0405060201 


0.99 


pn 


231 


l.lle-02 










MOS 


214 


7.48e-03 






0405060301 


0.44 


MOS 


107 


8.32e-03 




WR79 


5372 


0.72 


ACIS-S 


384 


5.04e-03 






6291 


0.23 


ACIS-S 


203 


4.57e-03 





Note - The binary phase for each observation was calculated using the value for To from Pourbaix et al. (2004; The ninth catalogue of 
spectroscopic binary orbits) except for WR46 (see notes to Table[JJl. The orbital 'phase' for WR145 is with respect of the first in time 
X-ray observation since no value for To is available in the literature. The moment To (binary phase 0.0) is the time of maximu m value 
of the radial velocity of the WR component in the binary system. The values for the orbital period in [Pourbaix et al ] [|2004r i arc 
identical to those in Table [JJ except for WR46 and WR145. For these two, we used more recent sources given in the note s to Table U 
The last co lumn gives the references where the c orresponding X-ray d ata h ave been already discu s sed in som e detail: (l)lBhatt et alj 
fcOld) ; (2) lFauchez et al.1 1I2OHII : (3) - iGosset et al] <2011aD ; (41 iHenault-Brunet et al] ll201ll l; (51 iNazel d2009h . 



absorption higher than that expected from the optical/UV 
observations. Since the interaction region in close CSW bi- 
naries should be located near the massive stars, it is realistic 
to attribute this extra X-ray absorption to the massive stel- 
lar winds of these objects. For sake of simplicity, we assumed 
that the extra absorption is due to the stellar wind of the 
WR star in the system because it is considerably more mas- 
sive than that of its companion. 



Thus, the models for our final global spectral fits con- 
sisted of the following components. An absorption compo- 
nent due to the interstellar matter towards the studied ob- 
ject (model wabs in XSPEC). The value of the neutral hy- 
drogen column density was kept fixed to that based on the 
optical/UV obser vations (e.g. A„ in Table [T]) through the 
iGorensteinl (119751 1 conversion (N H = 2.22 x 10 21 Av- cm" 2 ). 
Anticipating the discussion of the derived results, we note 
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that our conclusions from this study do not change if we use 
a mo re recent conversion |Vuong et al.ll2003l ; iGetman et al.l 
120051 : N H = 1.6 - 1.7 x 10 21 Ay cm" 2 ). Emission compo- 
nents are represented with discrete-temperature optically- 
thin plasma models (apec in XSPEC) with individual 'wind' 
absorption for each plasma component. The emission and 
'wind' absorption components had the same abundances. 
We recall that in the framework of the CSW picture, 
the discrete-temperature models represent qualitatively the 
te mperature str atification of the interaction region (see §5.2 
in IZhekovll2007h . We note that there are no objects in our 
sample of close CSW binaries that have elliptical orbits. This 
means that the amount of hot gas in the interaction region 
remains the same over an orbital period. So, for each ob- 
ject that had multi-phase observations, we adopted the same 
model normalization for the corresponding plasma compo- 
nent. This is equivalent to a non- variable emission measure 
with the orbital phase that is there is no variation in the 
intrinsic X-ray emission. We thus attribute a possible X-ray 
variability to the variable X-ray absorption with the orbital 
phase. 

One-temperature model was successful for fitting the 
X-ray spectrum of WR141. We found that two-temperature 
models gave adequate global fits for the X-ray emission from 
all other WN and WC stars in the sample but WR46. For 
the latter, we needed a three-temperature model. We note 
that the two-temperature one gave a statistically acceptable 
result but it was not able to match all the spectral details 
in the spectrum: e.g., it produced almost none emission for 
the Si XIII line complex at ~ 1.85 keV and variable silicon 
abundances could not fix this problem. It is our understand- 
ing that a more complex model was needed in this case due 
to the good photon statistics of the WR46 spectra (the best 
quality data in our sample). Also, we had no data on the in- 
terstellar absorption towards WR145, we thus fitted for the 
X-ray absorption in the global spectral fits for this object. 

Figure[T]and Table present the results from the global 
fits to the X-ray spectra of our sample of close WR+O 
binaries. We see that the quality of the fits is acceptable 
but we note that in many of the studied cases the for- 
mal quality of the fit (x 2 /dof) is good likely due to the 
low photon statistics. Nevertheless, we note that the fit re- 
sults firmly establish that a considerable amount of hot- 
temperature plasma (kT > 2 keV) is present in the studied 
objec ts. The 'canonical' WN an d WC abundances adopted 
here l|van der Hucht et al.lll986h seem to be an adequate ap- 
proximation to the actual abundances in the studied ob- 
jects. Although, a by-product from these fits in this respect 
is worth mentioning. Namely, as seen from Fig. [T] the SXV 
line complex at ~ 2.45 keV is not well matched by the cur- 
rent model. Our attempts to solve this problem by adding 
an additional temperature component were not successful. 
The only way to match the SXV 'bump' in the spectra was 
to adopt a variable sulfur abundance. Interestingly, the val- 
ues for the S abundance derived from all fits are not widely 
spread but 'cluster' around a factor of three increase with 
respect to the one in the 'canonical' WN and WC abun- 
dances. More specifically, the average sulfur abundance from 
all data sets is 2.97 ±0.33 (where the error is the error of the 
me an and not the varian c e) . W e note that the sulfur value in 
the Ivan der Hucht et al.l j 19861 ) WN and WC abundances is 
based on the cosmic mass fraction (see the footnote to their 



Table 1), thus, we get an indication from the current anal- 
ysis that the sulfur abundance in the 'canonical' WN and 
WC sets of abundances may need be increased by a factor 
of three. 

In the next section, we will discuss the global fit results 
in the framework of the CSW picture. 



5 DISCUSSION 

Although our sample of close WR+O binaries is limited, 
the adopted uniform way of analysing the X-ray emission 
from the studied objects provides useful pieces of informa- 
tion to check the global consistency of the colliding stel- 
lar wind picture in such objects. We note that our spec- 
tral analysis adopted simplified models: those of discrete- 
temperature optically-thin plasma (§ [4}. This means that 
the X-ray emitting plasma is in collisional ionization equi- 
librium (CIE), that is the non-equilibrium ionization (NEI) 
effects do not affect its emission. Also, the electron and ion 
temperatures are equal. We can use the corresponding char- 
acteristic plasma timescales to check validity of these as- 
sumptions in the case of close WR+O binaries. 

For this purpose, two di mensionless parameters have 
been introduced r eo = r. 



( Zhekov fc Skinned [2000 ) and 



NEI 



t/tnei (|Zhekovll2007h , where r is the timescale of 
the gasdynamics, r eq is the electron-ion temperature equal- 
ization time and tnei is the representative NEI timescale. 
We recall that the electron and ion temperatures are dif- 
ferent if T eq ^ 1 but their difference can be neglected if 
T eq 2> 1. Similarly, the NEI effects must be taken into ac- 
count if Tnei ^ 1 but can be neglected if Tnei S> 1. 

For a total mass of the WR+O systems of 30- 
40 Mq, we can use the Kepler's third law (a = 2.928 x 
10 11 P d 2/3 [M/M Q ] 1/3 cm; P d is the orbital period in days 
(Table (TJ; M is the total mass) to estimate the typical 
lengthscale for the gasdynamic problem (e.g., the binary 
separation). Then the stellar wind parameters from Ta ble [T] 
(M, Voo) along w i th eq .(l) in lZhekov fc Skinned (|2000h and 
eq.(l) in IZhekovl (|2007T ) show that, first, the plasma in the 
CSW shock in the close WR+O binaries considered here has 
equal electron and ion temperatures (F e? 3> 1). Second, the 
NEI effects are not important in these objects (Tnei 3> 1). 
All this makes us confident that the adopted spectral mod- 
els in this study were adequate to the physical conditions of 
the X-ray emitting plasma in close WR+O binaries. 

On the other hand, when fitting the total ('inte- 
grated') X-ray emission from temperature stratified plas- 
mas (e.g., present in the interaction region in CSW bi- 
naries) with discrete-temperature models, we can expect 
that the derived plasma temperatures will not be higher 
than the maximum one in the X-ray emitting region. We 
note that for a strong shock (with adiabatic index 7 = 
5/3) in helium-dominated plasma (mean molecular weight 
per particle fi = 4/3) the postshock temperature is kT 
= 2.608(Kh/1000 km s" 1 ) 2 keV, where V sh is the shock ve- 
locity, and the maximum shock velocity is equal to the stellar 
wind velocity (V B h = Voo) in the framework of CSW model. 
We thus see from Tables [1] and [3] that the deduced plasma 
temperatures in our sample of close WR+O binaries are 
qualitatively consistent with the adopted physical picture. 
Namely, all of them are lower than the maximum possible 
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Figure 1. The background-subtracted spectra overlaid with the best- fitting model. Each panel is marked by the name of the object and 
the binary phase in parentheses (see Tablcs[2]and [3). The spectra are re-binned to have a minimum of 10-20 counts per bin. In the panels 
with two different spectra in the same data set, the upper and lower (in red) curves are for the pn and MOS spectrum, respectively 



plasma temperature in the interaction region of these CSW 
objects. 

Also, in the CSW picture the hotter plasma is located 
near the line connecting the two stellar components in the 
binary system, that is in the denser part of the stellar winds. 
Thus, we can expect that the higher-temperature component 
of the discrete-temperature models will suffer higher 'wind' 



absorption. As seen from Table [3] this is in general found in 
the spectral fits to the X-ray emission from the close WR+O 
binaries in our sample (although there are some exceptions 
but likely due to the poor photon statistics). 

It then seems conclusive that the X-ray emission from 
the close WR+O binaries considered in this study is qual- 
itatively consistent with the physical picture where X-rays 
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Table 3. Global Spectral Model Results 



Name Phased X 2 / dof N H, ISM W Nh^ i (c) kTi^ 



EMi( £ ) 



N, 



,(<0 kT2 W 



EM 2 (e) Fx (/) 



WR46 



0.55 



231/255 2.10 



0.002 



+0.34 
0.002 



WR47 


0.62 
0.89 


244/295 


7.92 


0.53 
0.80 


WR139 


0.47 
0.55 
0.81 


362/335 


5.86 


0.23 
0.43 
0.10 


WR141 


0.02 


35/39 


8.20 


4.40 


WR145 


0.33 


180/278 


35.4 + || 


0.0 




0.00 




47 5 +5:5 


0.0 




0.09 




57 5+ 8 ' 4 
°'-°-6.9 


0.0 




0.18 




47 5+ 7 ' 5 


0.0 


WR148 


0.99 
0.44 


54/101 


5.16 


1.98 
1.03 


WR79 


0.72 
0.23 


40/49 


3.08 


0.09 
0.09 



+0.08 
-0.08 
+ 0.11 
-0.10 
+0.16 
-0.16 
+0.12 
-0.15 
+0.09 
-0.10 
+0.59 
-0.53 



+ 1.20 
-0.73 
+0.49 
-0.50 
+0.01 
-0.01 
+0.02 
-0.02 



0.15 



-0.01 
-0.01 



0.59 



+0.03 



0.58 



+0.04 
0.03 



3 40+ ' 43 

99+ ' 14 
u - 3a -0.12 



0.40 
1.27 



-0.18 
-0.11 

+0.48 
-0.20 



5.55 



+31.5 
0.27 



12.31^ 



+0.60 



3.66 



fi o 9 +0.75 
"■^-0.68 

16.5+ 3 j 



44.5 



+ 123. 
26.5 



0.09 



+0.02 



+0.09 
0.09 
+0.72 
0.44 
+0.21 
0.00 
+0.78 
0.66 
+0.63 
0.58 
+0.29 
0.26 
+0.19 
0.18 



0.35 
0.44 
0.00 
4.33 
6.11 
2.07 
1.00 

0.0 
0.0 
0.0 
0.0 

J - Da _1.69 

10 i+ 52 - 4 

lu - x -6.43 
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13.4 (27.3) 
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6.32 (39.7) 
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4.59 (39.7) 
5.10 (39.7) 
0.255 (4.14) 
0.249 (4.14) 
1.80 (33.6) 
1.70 (33.6) 



Note — Fit results from the two- temperature optically-thin plasma model with individual wind absorption for each component. A 
three-temperature plasma model was used only in the case of WR46: parameters of the third component are given in the second line for 

this object. The uncertainties are la errors from the fit. 
( a )The binary phase for each observation from Table [2] 
( b 'The interstellar absorption column density in units of 10 21 cm -2 . 
' c 'The helium-dominated wind absorption column density in units of 10 21 cm~ 2 . 
( d )The plasma temperature is in keV. 
( e )The emission measure (EM = J n e njj e dV) in units of 10 54 cm -3 at the distance for each object from Table[T] 
'^'The observed flux (0.5 - 10 keV) followed in parentheses by the unabsorbed value. The units are 10 — 13 erg cm -2 s — 1 . 



arise in hot temperature-stratified plasmas behind colliding 
stellar wind shocks. But, an interesting issue is worth dis- 
cussing as well, namely, the total X-ray energetics (luminos- 
ity) of CSWs in close binaries which is directly related to 
whether the CSW shocks are radiative, adiabatic etc. 

We can u se the dimens i onless parameter \ — t coo i/t 
introduced by IStevens et all l|l992l ) to estimate the impor- 
tance of radiative losses behind the CSW shocks in the 
WR+O binaries of our sample. We recall that r coo ; and r are 
the characteristic cooling time of the shocked plasma and the 
timescale of the gasdynamics, respectively. Then from the 
data in Table [T] and eq.(8) in IStevens et al.l (|l992l ). we see 
that the radiative cooling should have an important effect 
on the physics of CSWs in the close WR+O binaries studied 
here (x < 1)- This means that only numerical hydrodynamic 
modelling of the interaction region must be used to confront 
theory and observations. However, we can perform a qual- 
itative check on the X-ray energetics, as mentioned above, 
by considering two extremes: (a) highly radiative and (b) 
adiabatic CSW shocks. 

In the case of radiative shocks, we have an upper limit 
on the available energy (luminosity) that can be converted 
into X-ray emission: no more energy is emitted than the 
energy flux crossing the shock front per unit area. For the 
maximum X-ray luminosity, we have Lx = | J pV^ ind ^dS 
(p is the density of the wind in front of the shock; V w i n d,± 
is the wind velocity component perpendicular to the shock 
front; S is the shock surface). In the case of CSWs, the inte- 
gration is over the entire CSW 'cone' and although its shape 
varies from one object to another we can expect that a rela- 



tion Lx oc L W i n d might hold in general for highly radiative 
CSWs (L winc i = iMV^, is the stellar wind luminosity). 

On the other hand, in the case of adiabatic CSWs 
there exists a scaling law for the CSW X-ray luminosity 
with the mass- loss rate (M), wind veloci ty (Voo) and bi- 
nary separation (a): Lx o c M 2 V ^ 3 a _1 l|Luo et alj Il990l : 
iMvasnikov fc Zhekovl Il993l ). So, using the third Kepler's 
law we can write: Lx oc C'csw by introducing Ccsw = 
(Af/KT 5 M Q yr- 1 ) 2 (l/ O o/1000 km s- 1 )- 3 P d " 2/3 as a 'col- 
liding stellar wind' parameter. 

Figure [2] presents two plots: the X-ray vs. the wind lu- 
minosity and the X-ray luminosity vs. the CSW parameter 
for the objects in our WR+O binaries sample (the unab- 
sorbed fluxes from Table are used for calculating the X- 
ray luminosity). We see that there is no correlation for the 
former while the X-ray luminosity clearly correlates with 
the CSW parameter. Interestingly, the derived proportion- 
ality between log Lx and log Ccsw is 1.05 + 0.09 (the error 
is la error from the fit), thus, Lx oc C^'g^ 0,09 . All this is 
quite surprising since it means that the CSW shocks in close 
WR+O binaries are more likely adiabatic rather than being 
strongly radiative. We note again that our sample of studied 
objects is limited but nevertheless the results in Fig. [2] are 
very interesting and we will discuss them in some detail. 

It is worth noting another result from Fig.[2]that favours 
the case of adiabatic CSWs in the close WR+O binaries in 
our sample. Namely, the low efficiency of converting the wind 
luminosity into X-ray emission: Lx/L W i„d ~ 10~ 4 — 10 -5 . 
Of course, this ratio cannot be very close to unity because 
the CSW shock cone occupies only part of the 'sky' of the 
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Figure 2. Left panel: X-ray (Lx) vs. wind luminosity (L w i nc i = jMV^,) for the sample of close CSW binaries. Right panel: X-ray 

luminosity vs. the CSW parameter, C C SW = (Af/10 -5 Mq yr- 1 ) 2 (V O o/1000 km s _1 )~ 3 P^ 2/3 . The slope of the solid line is 1.05±0.09. 
The X-ray luminosity is calculated for the corresponding distance to each object (Table [TJ. The error bars on Lx are internal only and 
do not take into account those on distance to the objects. 



WR star (the shocked WR wind dominates the X-ray emis- 
sion from the interaction region). But, it seems unlikely to 
explain such a low value of Lx / L m i„d only by geometry ef- 
fects and some other explanation is needed. Moreover, by its 
very definition the term 'radiative shocks' means that high 
percentage of the energy influx is radiated away (mostly in 
X-rays for fast shocks), while the opposite is valid for 'adi- 
abatic shocks'. Figure [3] shows the result from our numeri- 
cal integration for the maximum possible X-ray luminosity, 
Lx = | J pV^ ind ± dS, in the case of highly radiative CSW 
shocks. In such a case, the interaction region 'collapses' and 
its shape coincides with that of the contact discontinuity be- 
tween the shocked WR and O star winds. We adopted the 
ICanto et al. | (|l996t ) solution for the shape of the CSW ' cone . 
We see that the theoretical Lx/L w i n d values are consider- 
ably higher than those observed. Thus, the relatively low 
values for the X-ray luminosity, with respect to the wind lu- 
minosity of the more massive star in the binary, do support 
the conclusion that the CSWs in our sample of close WR+O 
binaries are adiabatic. 

But, adiabatic CSW shocks are possible in this case only 
if the mass-loss rate of the stellar winds in the binaries is 
considerably smaller, e.g. at least by an order of magnitude, 
than the values listed in Table [T] (we note that in such a case 
the actual values of the CSW parameter will be smaller and 
the plot in the right panel of Fig. [2] will be shifted to the left 
in the x-axis; a detailed analysis in that case is needed: see 
discussion below) . In other words, our analysis of the X-ray 
emission from close WR+O binaries finds an indication of 
considerable clumping in the winds of these massive stars. 
It then seems plausible to propose that the stellar wind of a 
massive WR (O) star is likely a two-component flow: a rar- 
efied continuous component with low mass-loss rate and a 
'discrete' component consisting of numerous dense clumps. 
The first component forms the CSW interaction region, thus, 
it plays role for the X-ray emission from close WR+O bina- 
ries while the second one dominates their optical/UV spec- 
tra. We note that in the last two decades the physical picture 
of clumpy stellar winds in the hot massive stars gathered ob- 
servational evidence and has become a standard ingredient 
in the analysis o f the optical, UV and infrared spectra of 
these objects (see lPuls et alj|2008l for a recent review on the 



matter). Here from the analysis of X-ray spectra of close 
WR+O binaries, we find additional evidence that the stel- 
lar wind in massive stars is clumpy (highly inhomogeneous) 
on a lengthscale at least of the order of the binary separa- 
tion (approximately the size of the CSW region). It is worth 
noting though that we need more detailed knowledge, both 
observational and theoretical, in order to build a coherent 
physical picture of clumpy stellar winds in massive stars. 
For example, what is the origin of the clumps; what is the 
distribution of clumps by their size; what is the clump evo- 
lution with the distance from the star; do the homogeneous 
component of the wind and the clumps share the same bulk 
velocity; thus, what is the efficiency of the radiative force 
that drives both of them; do both components of the stellar 
wind reach terminal wind velocity in front of the interaction 
region; does radiative braking play an important role for this 
complex wind structure? 

An interesting issue in this respect is about the fate of 
the dense clumps when interacting with the CSW region of 
the rarefied homogeneous components of t he ste llar winds 
in close WR+O binaries. ICherepashchukl £l990) proposed 
a qualitative picture where the dense clumps pass freely 
through the interaction region and part of them are deceler- 
ated in the O-star photosphere where their X-rays a re trans- 
forme d into optical emission. On the other hand, iPittardl 
(2007) carried hydrodynamical simulations of a CSW bi- 
nary with clumpy winds that showed efficient destruction of 
clumps while interacting with the CSW shocks. As a result, 
the average density in the interaction region increases and is 
similar to the case of smooth winds with the higher mass-loss 
rate. We have to keep in mind that these simulations were 
suitable for wide CSW binaries and specificall y for the stel- 
lar wind and binary parameters of WR140 (see |Pittardll2007l 
for details) and the reality might be quite different in close 
CSW binaries. Namely, the clumps are likely smaller closer 
to the base of the stellar wind, they may have much larger 
density contrast with respect to the smooth component etc. 
Also, if the clumps are destroyed in the CSW region in close 
WR+O binaries, then these objects would have been much 
more luminous in X-rays than deduced from observations. 
Thus, it seems plausible to assume that the dense clumps of 
the stellar winds 'survive' while crossing the CSW region in 
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Figure 3. Left panel: the ratio of the maximum possible X-ray luminosity (Lx) to the wind luminosity (L w i n d = ^MV^) for 
different values of the half-opening angle (/?) of the CSW 'cone' in the case of highly radiative CSW shocks in binary system. Right 
panel: dependence of the half-opening angle (/3) on the ratio of the ram pressure of the stellar winds in the binary system (A = 
MwrVwr/MoVo, where Mwr,V\yr, Mq,Vo are the wind parameters, mass loss rate and wind velocity, of the WR and O star, 
respectively). 



close WR+O binaries. It is difficult to describe what could 
be the possible observational evidences from the interaction 
of these dense clumps with the CSW region since they de- 
pend on such 'hard-to-guess' details as those mentioned at 
the end of the previous paragraph. 

We thus believe that future X-ray observations with 
much better photon statistics might be very helpful in this 
respect. Results from such observations must be considered 
in conjunction with those from the optical/UV spectral do- 
main to build a self-consistent physical picture of the stellar 
winds in massive stars. Such a global analysis may reveal 
that the smooth component has very little or no contribu- 
tion to the optical/UV emission of a massive star and X-ray 
observations could be the only tool to reveal its presence 
and physical properties. It should also take into account 
that massive O s tars are X-ray sourc es themselves (for a 
recent review see iGiidel fc Naie1l2009l ). We just note that 
their emission is rather 'soft' with plasma temperatures be- 
low 1 keV fe.g.. lWoidowski fc SchuQl2005l ; IZhekov fc Pallal 
120071; see also Sections 4.1.3 and 4.3 in the review paper 
of lGiidel fc Naze 1120091 ). opposite to what is found in close 
WR+O binaries (e.g., Table [3]). On the other hand, the is- 
sue of the intrinsic contribution of the WR star to the total 
X-ray emission from the binary might still bear a lot of un- 
certainties. 

Pointed observations with modern X-ray observatories 
(Chandra; XMM-Newton) of a few p resumably single WC 
stars resulted only in non-detections (|Oskinova et al.ll2003l ; 
ISkirmer et al.| [2006). Could it be that all single WC stars are 
X-ray quite? In such a case, there will be no contribution 
from a WC star to the total X-ray emission from a WR+O 
binary system. 

In contrast to this, the lSkinner et ail (|201dh analysis of 
the X-ray spectra of a small sample of presumably single WN 
stars showed that these objects emit X-rays that arise from 
an admixture of cool (kT < 1 keV) and hot (kT > 2 keV) 
plasma. Presence of hot plasma (kT > 2 keV) makes them 
similar to the CSW binaries studied here and distinct from 
the single O stars. So, could it be that these WN stars are 
not single but binaries instead? However, there are some dif- 
ferences in the X-ray characteristics between the single WN 



stars and the close WR+O binaries in our sample. Namely, 
the latter are X-ray more luminous with a minimum value 
of log Lx ~ 32.5 erg s^ 1 (Fig. [2| while most of the X-ray de- 
tected presumably single WN stars have lumin osities below 
this figure (see Fig. 10 in lSkinner et al.l 120101 1. Also, there 
is indication that the Lx oc L min d relation holds for sin- 
gle WNs while this is not the case for c lose WR+O binaries 
(compare Fig. 10 in lSkinner et alj|2010l with the left panel in 
Fig. [5]here). Since the X-ray production mechanism in single 
WN stars has not been identified yet, one possibility is that 
those of them detected in X-rays are in fact WR binaries 
with a nor mal (non-degenerat e) companion (see discussion 
in § 4.4 in Skinner et al.l 2010) which may also explain the 
Lx oc L W i n d trend for them (see § 7.2 in IZhekov fc Parkl 
l2010bl ) . And we emphasize again that more X-ray data with 
good quality are needed and the global analysis as that men- 
tioned above will help us better understand the physical pic- 
ture in the stellar winds of hot massive stars. 

Finally, the object with the shortest period in our sam- 
ple, WR46, deserves a few more comments. This object has a 
very complex v ariability pattern in the optical, UV and may 
be X-rays (e.g.. iMarchenko et al.ll2000l ; iGosset et ai"1l2011al ; 
i Henault-Brunet et alj|201ll ; but see also Appendix lAl here) . 
Icosset et al.l l|2011al ) prop osed that its hard X-ray emis sion 
may arise in CSW shocks. IHenault-Brunet et al.l (|201lh ar- 
gued for non-radial pulsations as the most likely scenario 
to explain its characteristics and they even drew analogy 
between WR46 and C, Pup (a massive presumably single O 
star) in this respect. Based on the results presented here 
(e.g., high plasma temperatures) and since the X-ray vari- 
ability pattern is not well established (see Appendix 0, we 
believe that the physical picture of CSWs in a short-period 
WR+O binary is a more likely explanation for the X-ray 
properties of WR46. 



6 CONCLUSIONS 

Using data from the Chandra and XMM-Newton public 
archives, we have analysed the X-ray emission from a small 
sample of close WR+O binaries. In such objects, X-rays 
likely originate in colliding stellar wind shocks driven by the 
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massive winds of the binary components. The main results 
and conclusions from our analysis are as follows. 

(1) Global spectral fits show that two-temperature 
plasma is needed to match the X-ray emission from these 
objects as the hot component (kT > 2 keV) is an important 
ingredient of the spectral models. 

(2) In general, CSW shocks in close binaries are ex- 
pected to be radiative due to the high density of the plasma 
in the interaction region and the X-ray emission is domi- 
nated by that of the shocked WR wind. Thus, a correla- 
tion between the X-ray luminosity from these objects and 
the mechanical luminosity of the WR wind should exist: 
Lx oc L win d- Interestingly, we do not find such a correlation 
for the objects studied here. 

(3) Our analysis shows that a correlation between the 
X-ray luminosity and the so called CSW parameter (see §[S]) 
does hold. This means that CSWs in close WR+O binaries 
must be adiabatic and this is possible only if the mass-loss 
rates of the hot stars in these objects are at least one order 
of magnitude smaller than the values currently accepted. 

(4) The most likely explanation for the X-ray proper- 
ties of close WR+O binaries could be that their winds are 
two-component flows. The more massive component (dense 
clumps) play role for the optical/UV emission from these 
objects. However, the smooth rarefied component is a key 
factor for their X-ray emission. We believe that global anal- 
ysis (modelling) of optical, UV, X-ray emission from close 
WR+O binaries will help us build a self-consistent physi- 
cal picture of the stellar winds and the close circumstellar 
environment in these objects. 

(5) To further check the results presented here, simi- 
lar analysis should be applied to close O+O binaries since 
they are evo lutionary progeni t ors o f the WR+O binaries. 
The work bv lDe Becker et al.l (|2004T l is an example for the 
relevance of such a proposition. Their study of the X-ray 
emission from HD159176 (an 0+0 system with a 3.367-day 
period) showed that the CSW model overestimates the ob- 
served X-ray luminosity for the standard wind parameters. 
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APPENDIX A: DETAILS ON THE X-RAY 
DATA FOR WR46 AND WR79 

Al WR46 



One of the most interesting 
is its variability established m 



characteristics of WR46 
_ the optical/UV (e.g., 
iMarchenko et al"]|2000l and the references therein). Since 
XMM-Newton is capable of providing optical/UV data si- 
multaneously with the X-ray data, the corresponding WR46 
observations allow to search for X-ray variabil ity that may 
correla te wi th the one found in the optica l/UV. iGosset et all 
(|2011ah and lHenault-Brunet et all (|2011 ) reported that such 
a correlation does exist (we note that Gosset et alj|2011al 
suggested that the X-ray variability is restricted only to the 
0.2-0.5 keV energy range). 

Adopting a recent version of the SAS software ( § [3J , 
we extracted background-subtracted X-ray light curves (LC) 
from the reprocessed EPIC data. We also built UV LCs in 
the two filters (UVM2, UVW2) from the reprocessed data 
of the Optical Monitor used in these observations. For con- 
sistency with the UV data, the X-ray LCs are binned in 
2000-second time intervals. Figure I All presents the corre- 
sponding X-ray and UV LCs of WR46. We see that the 
variability pattern in WR46 is quite complex. Namely, the 
expected 0.329-day period is evident only in the UV filter 
UVW2. This is not the case for the other UV filter (UVM2) 
although these data are taken one after the other (see the 
left panel in the first row in Fig. I All) . Such a behaviour is 
quite strange (or interesting) since the throughput curves 
for these filters overlap considerably (see Fig. 88, § 3.5.5.1 



in the XMM-Newton Users HandboolQ). Similarly, the X- 
ray variability (if any) does not seem to be 'identical' in the 
EPIC detectors of a different kind: pn and MOS. 

To quantify the variability pattern in WR46, we per- 
formed the following exercise that is indicative of whether 
the corresponding emission is variable. First, we fitted each 
LC with a constant (adopting \ 2 fitting) . Second, excluding 
the two 'extreme' values, i.e. the maximum and minimum 
value, in each LC, we fitted the LC with a constant again. 
The result as a formal goodness of the fit is as follows. For 
the complete LC, the goodnes of the fit for a constant emis- 
sion is 0.0002 (UVM2), 0.0 (UMW2), 0.49 (X-ray, pn) and 
0.99 (X-ray, MOS). In the case with the two extreme val- 
ues excluded, the goodness of the fit is 0.24 (UVM2), 0.0 
(UMW2), 0.87 (X-ray, pn) and 1.0 (X-ray, MOS). We note 
that excluding one or two data points from a sequence of 
measurements is not expected to change the global trend in 
a variability pattern of a physical quantity. Thus, we feel it is 
safe to conclude that only the UV emission of WR46 in the 
UVW2 filter of the optical monitor on-board XMM-Newton 
does show clear sign for not being constant. We believe that 
more X-ray and UV data, taken simultaneously, are needed 
to firmly establish whether the X-ray emission from WR46 
is variable, and if it is, then on what timescale, on what lu- 
minosity scale and whether it correlates with the variability 
detected in other spectral domains (e.g., optical, UV). 



A2 WR79 

Th e XMM-New t on obse rvations of WR79 were analysed 
by IGosset et alj (|2011bh . But thanks to the Chandra su- 
perb spatial resolution, the ACIS-S images of WR79 re- 
veal that there are two near-by sources within 6" from 
WR79 (Fig. 1X2]) . We denote them as SI (the northern 
source) and S2 (the southern source). The SI coordinates 
are q 20 oo = 16 h 54 m 19!41, <5 2 ooo = -41°49'08'.'15 and WR79 
is the only object in SIMBAD within 5" from SI (radial dis- 
tance of 4'.'7). Source SI is thus unidentified yet. The S2 co- 
ordinates are q 2 ooo = 16 ,l 54 m 19 , :44, <5 2 ooo = -41°49'16'.'71 
(radial distance of 5'.' 9 to WR79) and the binary system 
CCDM J16543-4149D is the only object in SIMBAD within 
2" from S2 (radial distance of l'.'9). Source S2 is thus uniden- 
tified yet. 

It is important to note that it is possible to extract 
a 'net' X-ray spectrum of WR79 from the Chandra data 
while only a total spectrum for all three sources, WR79, 

51 and S2, can be obtained from the XMM-Newton data 
(Figs. rAJ and \K3&- We see from Fig. \M\ that SI is defi- 
nitely a variable source and Figure IA3I illustrates that the 
X-ray emission from SI and S2 may alter appreciably the 
X-ray emission from WR79 when the latter is observed with 
X-ray telescope having not very high spatial resolution. We 
note that the total number of counts from sources SI and 

52 amounts to 37-47% of that from WR79 as directly mea- 
sured in the ACIS-S data of the two Chandra observations. 
Also, the X-ray emission of SI and S2 is much softer com- 
pared to that of WR79 and if the total spectrum of these 
three sources is assigned to WR79 alone, then from analy- 
sis of such a 'combined' spectrum we may draw conclusions 



see Documents & Manuals on |http://xmm.esac.esa.int/ 
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Figure Al. XMM-Newton light curves (LC) of WR46. First row: LCs as a function of observing time in two UV filters of the Optical 
Monitor (UVM2; UVW2) and in the (0.2 - 10 keV) energy range with the EPIC detectors (pn; MOS denotes the total emission from the 
two MPS detectors). M iddle and Bottom rows: LCs folded with the orbital period; the values for period and zero-phase time are from 
iMarchenko et al.l feOOOf) . Note that for graphical convenience the data points are plotted twice: the data for phase values (1.0-2.0) are the 
same as for phase values (0.0-1.0). Circles mark the maximum and minimum values in each LC. The X-ray LCs are background-subtracted. 
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Figure A2. WR79 images: Chandra ACIS-S images for the two observations used in this study with ObsID 5372 (Obs 1) and 6291 
(Obs 2), and an XMM-Newton MOS1 image (ObsID 0109490401). All images are centered at the optical position of WR79 (SIMBAD) 
marked by a circle with crosshairs. Two X-ray sources are seen in close distance from WR79: marked by a cirle and denoted SI and S2. 
A square marks the optical position of CCDM J16543-4149D (a binary system; SIMBAD). The green circle has a diameter of 15". 
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about the X-ray characteristics (e.g., plasma temperature, 
X-ray absorptions, variability) that will not be quite correct. 
From all this, we can safely conclude that for the moment 
Chandra observations are the only ones that can provide us 
with valuable information on the X-ray emission from the 
WR+O binary WR79. 

This paper has been typeset from a TpX/ WF^K file prepared 
by the author. 
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Figure A3. The background-subtracted spectra of the X-ray sources from Fig. IA2I The WR79 spectra extracted from the Chandra 
ACIS-S data are shown in the first two panels of the upper row as those from the large extraction region having a diameter of 15 arcsec 
(shown in green in Fig. I A2ft are drawn in red while the spectra only of WR79 itself are given in black (they are the ones used in the 
current study, see Fig.[TJ. For comparison, the archive XMM-Newton MOS1 spectrum of WR79 is shown in the last panel of the upper 
row. The X-ray spectra of sources SI and S2 (Fig. I A2t are shown in the lower row. No spectrum of SI from Obs 2 is presented since 
there are only 2 counts in it. All the spectra are re-binned to have a minimum 10 cts per bin. 



